Impact of Optical Modes on the Pairing Potential in Bilayer Cuprates 



LA. Larionov* and M.V. Eremin 

Magnetic Radiospectroscopy Laboratory, Department of Physics, Kazan State University, 420008 Kazan, Russia 

Superconducting transition temperature is calculated for differently doped bilayer cuprates. Su- 
perexchange is assumed to be the dominating mechanism of high-temperature superconductivity, 
but the contribution from the phonon potential is not negligible, which qualitatively explains the 
observed weak isotopic effect. The calculated value 2A ma x /ksTc — 4.5 is close to the experiment 
in the case of optimum doping. 
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I. INTRODUCTION 

The High Temperature Superconductivity (HTSC) was discovered 11 years ago. None of the proposed HTSC 
mechanisms has been commonly accepted yet (see [1]). Unlike the case in normal low-temperature superconductors, 
the energy gap in HTSC materials is different at different points of the Brillouin zone. This dependence is important, 
for the understanding of the HTSC mechanism. 

In [2,3] we found solutions of the Bardeen-Cooper-Schrieffer (BCS) equations taking copper spin superexchange 
into consideration. This interaction was found to lead to a reasonable critical temperature Tq and d-symmetry of the 
energy gap in superconducting cuprates. The conclusion about the d-symmetry of the order parameter was confirmed 
by recent studies of photoemission spectra [4]. In the present work we additionally consider the electron - phonon 
interaction with optical buckling A\ g , B\ g and breathing A g , E g modes to explain the weak isotope effect in these 
compounds. 

II. POTENTIAL OF CURRENT CARRIERS PAIRING THROUGH OPTICAL BUCKLING AND 

BREATHING MODES 

In the compound YBa2Cu307-a; holes interact most strongly with an electric field perpendicular to the CuC>2 plane. 
This field is mainly induced by triply charged yttrium ions [5]. The operator of binding energy with oscillations 
!h ' perependicular to the CuC>2 plane has the form 

a ; 

H h - ph = e ]T {E x u x (an + ^) p<Z + E y u y (an + ^) p™}, (1) 
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where p™ and p™ are the Hubbard operators of oxygen holes, u x (n) and u y (n) are the displacements vectors of 0(2) 
and 0(3) positions in a unit cell with the number n, x and y are unit vectors of the axes a and b respectively, a is the 
lattice constant, and E x = 1.2 • 10 8 V-cm -1 , E y = 1.5 • 10 s V-cm -1 are the components of electric field along the the 
axis c, calculated in [6]. For simplicity, we put E = E x = E y = 1.35 • 10 8 V-cm -1 . Then, calculating the commutator 



X 



and using the expression 
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d , H h . ph ] = £ U«(q)*^ q (fo q + &t q ), 



where \I/j c ' p is the quasiparticle operator of singlet-correlated oxygen holes, and (6 q are the phonon creation and 
annihilation operators, one finds 
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where u Q = wa 1s =440 cm -1 and wb 1s =340 cm -1 [7]. The plus and minus signs correspond to the modes Ai g and 
Bi g respectively, e is the electron charge, and m is the mass of the unit cell. Since the difference between wa 1s and 
u>B lg is not so important, they are put equal to ujq =400 cm -1 . 

Using the Frohlich procedure, the potential of current carriers interating through the phonon field of buckling modes 
Ai s and Bi 9 is written as 



G(k' - k) = 2G f + h[^ (K ~k x ) a + cos (k^ y -k y ) a] h ^ 
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where Go = 35 meV. 

Assuming that the constant of interaction with the modes A g and E g (and their frequencies) are identical, the 
contribution to the pairing potential is written as 



B(k' - k) = 2B 2 
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According to the estimate of [8] , the constant B is about 60 meV. The frequencies wa s and we s , taken from [9] , are 
about 480 cm -1 . 



III. EQUATION FOR THE GAP 



The BCS equation for the energy gap, taking into account the interaction via the phonon field and the superexchangc 
between copper spins, has the form 



Ak ' = E { p2 t G ( k ' - k ) + B ( k ' - k )i ( hujD - i £k ' - ek i) - Jdd{ ^ - k )} (* k ' pd *lf ) , 
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where J (q) is the Fourier transform of the superexchange parameter J = 57 meV, 

J (q) = 2 J (cos q x a + cos q y a) , 
oj d = 500 cm -1 is the Debye frequency, and P = 5(1 + x). The correlation function is given by 
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where £4. = \J (e k — l^) 2 + The chemical potential \x and the hole concentration x for two copper sites in the 
bilayer (assuming that the antibonding band is ampty) are related as 



At x = 0.33 the chemical potential is placed on 10 meV below the saddle singularity peak in the density of states. 
The dispersion e k is chosen as 



e k = P [2ti (cos k x a + cos k y a) + At% cos k x a cos k y a + 2t% (cos 2k x a + cos 2k y a)] 



(6) 



where t\, ti, and £3 are the effective hopping integrals. All the calculations are carried out at t\ — 70 meV, ti = 0, 
and £3 = 5 meV. 
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FIG. 1. Temperature dependence of the superconducting 
transition temperature in YBa2Cu307-i/ with doping: solid 
line (calculations) and the experimental points are sketched 
by open circles. 

Figure 1 shows the dependence of the superconducting 
transition temperature of YBa 2 Cu3C , 7_, / with doping. 
It is evident that the calculations qualitatively conform 
with the experimental data (normalized by the expression 
Tc/T c , ma x = 1 - 82.6(a; - 0.32) 2 [10]). However, their 
noncoincidence remains unexplained. Numerical solution 
of the energy gap equation yields a <i-type symmetry of 
the order parameter Ak = A (cos k x a — cos k y a) and a 
nonstandard value of 2A max /kBTc « 4.5 with the su- 
perconducting transition temperature Tc,max ~ 100 K. 
Note that solutions of (5) with J = lead to the s-type 
symmetry of the gap, Ak = Ao(cosfc x a + cosk y a). The 
conclusion of [11] that the cZ-symmetry of the gap arises 
from pairing only through optical buckling modes seems 
to be groundless. 
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IV. CONCLUSION 



Self-consistent solutions of the BCS gap equation are 
found in the class of short-range pairing potentials for 
various chemical potentials. When the Fermi level ep is 
near the bottom(top) of the band, the solutions corre- 
spond to the s-type pairing, while for ep in the center 
of the band the solution are related to the d-type. The 
short-range potentials considered are (i) superexchange 
interaction, (ii) interaction of current carriers via opti- 
cal oscillations, and (iii) breathing and buckling modes 
of oxygen atoms in CuC>2 planes. The superconducting 
transition temperature is calculated for various oxygen 
indices x. 

The dominant HTSC mechanism is a superexchange. 
However, the contribution of the phonon pairing is not 
negligible, that qualitatively explains the observed weak 
isotope effect. The calculated value 2A max /kBTc « 4.5 
is close to the experiment. 
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